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We are presenting the first ab initio structure investigation of the looseiy bound lx Be nucleus, 
together with a study of the lighter isotope 9 Be. The nuclear structure of these isotopes is par- 
ticularly interesting due to the appearance of a parity-inverted ground state in 11 Be. Our study 
is performed in the framework of the ab initio no-core shell model. Results obtained using four 
different, high-precision two-nucleon interactions, in model spaces up to 9ftf2, are shown. For both 
nuclei, and all potentials, we reach convergence in the level ordering of positive- and negative-parity 
spectra separately. Concerning their relative position, the positive-parity states are always too high 
in excitation energy, but a fast drop with respect to the negative-parity spectrum is observed when 
the model space is increased. This behavior is most dramatic for 11 Be. In the largest model space 
we were able to reach, the l/2 + level has dropped down to become either the first or the second 
excited state, depending on which interaction we use. We also observe a contrasting behavior in the 
convergence patterns for different two-nucleon potentials, and argue that a three-nucleon interaction 
is needed to explain the parity inversion. Furthermore, large-basis calculations of 13 C and n B are 
performed. This allows us to study the systematics of the position of the first unnatural-parity 
state in the N — 7 isotone and the A = 11 isobar. The 11 B run in the 9hQ model space involves a 
matrix with dimension exceeding 1.1 x 10 9 , and is our largest calculation so far. We present results 
on binding energies, excitation spectra, level configurations, radii, electromagnetic observables, and 
10 Be + n overlap functions. 

PACS numbers: 21.60.Cs, 21.45.+V, 21.30.-x, 21.30.Fe, 27.20.+n 



I. INTRODUCTION 

Studies of how nuclear structure evolves when varying 
the N/Z ratio are important in order to improve our fun- 
damental understanding of nuclear forces. For this rea- 
son, research on light neutron-rich nuclei has attracted an 
increasing amount of theoretical and experimental effort 
since the advent of radioactive nuclear beams. The appli- 
cation of a standard mean-field picture to describe these 
few-body systems is questionable, and it is not surprising 
that substantial deviations from regular shell structure 
has been observed. The A = 11 isobar is of particular in- 
terest in this respect since it exhibits some anomalous fea- 
tures that are not easily explained in a simple shell-model 
framework. Most importantly, the parity- inverted l/2 + 
ground state of 11 Be was noticed by Talmi and Unna Q 
already in the early 1960s, and it still remains one of the 
best examples of the disappearance of the N = 8 magic 
number. 

Many theoretical studies of odd- A beryllium isotopes 
have already been performed using various models. A 
thourough review of the structure of unstable light nu- 
clei in terms of the shell model can be found in Ref. Q. 
Of particular interest is the study on unnatural-parity 
states of the A = 11 isobar by Teeters and Kurath [3j us- 
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ing a lhfl model space and the Millener-Kurath interac- 
tion with modified Os and sd single-particle energies. The 
halo structure of the 11 Be ground state was reproduced 
with the variational shell model by Otsuka et al |4| . They 
used Skyrme interactions and constructed multi-nucleon 
wave functions from a variational single-particle basis in 
a (0-l)Ml model space. Alternatively, the loosely-bound 
nature of the valence neutron in 11 Be can be treated ex- 
plicitly in a 10 Be + n picture with a Woods-Saxon po- 
tential, see e.g. Refs. y, @- Using a coupled-channels 
treatment, the authors of these papers found a significant 
overlap with excited-core states. Possible explanations 
for the parity inversion of the 11 Be ground state has also 
been investigated using the AMD+HF model 0, which 
is a combination of anti-symmetrized molecular dynam- 
ics with the concept of single-particle motion. An ex- 
tended version of the AMD framework was later used to 
study excited states of 11 Bc, and the existence of three 
negative-parity rotational bands was proposed @ . 

There are also several calculations involving different 
cluster models. In particular, a-clustcring has been con- 
sidered to play an important role in these systems. With 
this assumption as a starting point, K. Arai et al used 
an a + a + n model and obtained the ground state of 9 Be 
using the stochastical variational method, while several 
particle-unbound excited states were studied simultane- 
ously with the complex scaling method A similar 
a + a + Xn description was employed by P. Descouve- 
mont jioj in his study of possible rotational bands in 



2 



9_11 Be using the Generator Coordinate Method. His 
conclusion, however, was that the degree of a-clustering 
decreased with increasing mass, and consequently his 
model was not able to reproduce some of the anomalous 
properties of 11 Be. 

Unfortunately, due to the complexity of the problem, 
there has been no genuine ab initio investigation of n Be 
starting from realistic inter-nucleon interactions. There 
is no doubt that the cluster and potential models are very 
successful, and can provide reasonable explanations for 
many observations. Still, one has to remember that they 
rely upon the fundamental approximation that the total 
wave function can be separated into cluster components. 
Furthermore, the effective interactions used in all mod- 
els must be fitted to some observables for each individual 
case. On the contrary, a truly microscopic theory such as 
the Green's Function Monte Carlo (GFMC) method fill . 
or the ab initio no-core shell model (NCSM) [H |l3|, 
starts from a realistic inter-nucleon potential and solves 
the A-body problem, producing an antisymmetrized total 
wave function. It is a true challenge of our understand- 
ing of atomic nuclei to investigate nuclear many-body 
systems (A > 4) using such ab initio approaches. 

This paper represents an effort to fill this gap. Our 
study is performed in the framework of the ab initio 
NCSM, in which the A-body Schrodinger equation is 
solved using a large Slater determinant harmonic oscilla- 
tor (HO) basis. However, it is well known that the HO 
basis functions have incorrect asymptotics which might 
be a problem when trying to describe loosely bound sys- 
tems. Therefore it is desirable to include as many terms 
as possible in the expansion of the total wave function. 
By restricting our study to two-nucleon (WW) interac- 
tions, even though the NCSM allows for the inclusion of 
three-body forces, we are able to maximize the model 
space and to better observe the convergence of our re- 
sults. 

In Sec. |H| the framework for the NCSM will be briefly 
outlined, and the four different high-precision WW inter- 
actions that are used in this work will be introduced. 
Sec. IIIII is devoted to a presentation and discussion of 
our complete set of results for 9 ' n Be, with a particular 
focus on the position of the first unnatural-parity state. 
Concluding remarks are presented in Sec. II VI 



II. AB INITIO NO-CORE SHELL MODEL 

Applying the ab initio NCSM is a multi-step process. 
The first step is to derive the effective interaction from 
the underlying inter-nucleon forces, and to transform 
it from relative coordinates into the single-particle M- 
scheme basis. The second step is to evaluate and di- 
agonalize the effective Hamiltonian in an A-nucleon (Z 
protons and W neutrons) Slater determinant HO basis 
that spans the complete W max ftf2 model space. Finally, 
we can use the resulting wave functions for further pro- 
cessing. This section contains a short discussion on each 



of these steps. We stress that an important strength of 
the method is the possibility to include virtually any type 
of inter-nucleon potential. The four high-precision AW 
interactions that have been used in this study will be in- 
troduced in Sec. Ill CI A more detailed description of the 
NCSM approach, as it is implemented in this study, can 
be found in, e.g., Ref. [Tflj . 



A. Hamiltonian and effective interactions 

The goal is to solve the A-body Schrodinger equation 
with an intrinsic Hamiltonian of the form 
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where m is the nucleon mass and Vnnaj is the WW inter- 
action including both strong and electromagnetic com- 
ponents. As mentioned earlier, we will not use three- 
body forces in this study since we strive to maximize 
the size of the model space. By adding a center-of- 
mass (CM) HO Hamiltonian Pg M = T CM + Ug M (where 
Ug M = Amrt 2 R 2 /2, R = E^r/A), we facilitate the 
use of a convenient HO basis. The modified Hamiltonian 
can be separated into one- and two-body terms 
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The next step is to divide the A-nucleon infinite HO ba- 
sis space into an active, finite model space (P) and an 
excluded space (Q = 1 — P). The model space consists of 
all configurations with < W max M7 excitations above the 
unperturbed ground state. In this approach there is no 
closed-shell core; meaning that all nucleons are active. 

Since we solve the many-body problem in a finite model 
space, the realistic WW interaction will yield pathologi- 
cal results because of the short-range repulsion. Conse- 
quently, we employ effective interaction theory. In the 
ab initio NCSM approach, the model-space dependent 
effective interaction is constructed by performing a uni- 
tary transformation of the Hamiltonian J2J), e~ s H^e s , 
such that the model space and the excluded space are 
decoupled Qe~ s H^e s P = 0. This procedure has been 
discussed by Lee and Suzuki 0,0], and yields a Hermi- 
tian effective interaction P ff = Pe^ s H^e s P which acts 
in the model space and reproduces exactly a subset of the 
cigenspectrum to the full-space Hamiltonian. In general, 
this effective Hamiltonian will be an A-body operator 
which is essentially as difficult to construct as to solve the 
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full yl-body problem. In this study we approximate the 
effective Hamiltonian at the two-body cluster level. The 
basic idea is to derive it from high-precision solutions to 
the two-body problem with 7i 2 = h\ + h 2 + V^ A , where 
the two-body term is the same as in Eq. ©. The final 
result will be a two-body effective interaction V^'^j. See 
Ref. HUTU for details. 

We note that our approximated effective interaction 
will depend on the nucleon number A, the HO frequency 
fl, and the size of the model space which is defined by 
^Vniax- However, by construction, it will approach the 
starting bare interaction V i 
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Consequently, the dependence on ft will decrease with 
increasing model space, and the NCSM results will con- 
verge to the exact solution. A very important feature 
of the NCSM is the fact that the effective interaction 
is translationally invariant so that the solutions can be 
factorized into a CM component times a wave function 
corresponding to the internal degrees of freedom. Due to 
this property it is straightforward to remove CM effects 
completely from all observables. 



B. Solution of the many-body Schrodinger 
equation 

Once the effective interaction has been derived, we can 
construct the effective Hamiltonian in the A-body space. 
In this process we subtract the CM Hamiltonian Hq M , 
and add the Lawson projection term 0(Hcm — §M^) to 
shift eigenstates with excited CM motion up to high en- 
ergies. States with 05 CM motion are not affected by 
this term and, consequently, their eigenenergies will be 
independent of the particular choice of f3. We are now 
left with a Hamiltonian of the form 
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The computational problem of obtaining the many- 
body eigenvalues is non-trivial due to the very large 
matrix dimensions involved. The largest model space 
that we encountered in this study was the n B 9M7 
(positive parity) space, for which the dimension exceeds 
dp = 1.1 x 10 9 . For n Be( 9 Be), the 9 Ml space gives 
dp = 7.1 x 10 8 (2.0 x 10 8 ). To solve this problem we have 
used a specialized version of the shell model code AN- 
toine mm^ recently adapted to the NCSM This 
code works in the M scheme for basis states, and uses the 
Lanczos algorithm for diagonalization. The number of it- 
erations needed to converge the first eigenstates is signif- 
icantly reduced by the implementation of a sophisticated 
strategy for selecting the pivot vectors. This feature of 
the code is absolutely crucial when using it to perform 
calculations in very large model spaces. 



Furthermore, the code takes advantage of the fact that 
the dimensions of the neutron and proton spaces are 
small with respect to the full dimension. Therefore, be- 
fore the diagonalization, all the matrix elements involv- 
ing one- and two-body operators acting in a single sub- 
space (proton or neutron) are calculated and stored. As 
a consequence, during the Lanczos procedure, all non- 
zero proton-proton and neutron-neutron matrix elements 
can be generated with a simple loop. Furthermore, the 
proton-neutron matrix elements are obtained with three 
integer additions 01- However, for no-core calculations 
(in which all nucleons are active) the number of shells 
and, consequently, the number of matrix elements that 
are precalculated, becomes very large. Consider, e.g., 
the n B calculation in the 9M7 space. The full dimension 
is dp = 1.1 x 10 9 , while the number of active shells is 
66 and the total number of neutron plus proton Slater 
determinants is N(n) + N(p) = 1.0 x 10 7 . This corre- 
sponds to 80 Gb of precalculated and stored matrix ele- 
ments. In contrast, consider a shell-model calculation of 
57 Ni in the full fp space. The total dimension is larger, 
dp = 1.4 x 10 9 , but there are only four active shells which 
gives N(n) + N(j>) = 1.8 x 10 6 and requires merely 1 Gb 
of precalculated data. 

A recent development of the NCSM is the ability to 
further process the wave functions, resulting from the 
shell-model calculation, to obtain translationally invari- 
ant densities [2(| and cluster form factors [2lJ . The latter 
can be used to obtain spectroscopic factors, but can also 
serve as a starting point for an ah initio description of 
low-energy nuclear reactions. We have employed these 
new capabilities to study the overlap of 11 Be with differ- 
ent 10 Be + n channels. 



C. Realistic NN interactions 

Four different, high-precision NN interactions have 
been used in this study. These are: the Argonne V8' 
(AV8') jl E2, the CD-Bonn 2000 (CDB2k) the 
N 3 LO [2|, and the INOY [H m potentials. We can 
divide these interactions into three different types: 

1. Local in coordinate space: The AV8' interaction is 
an isospin-invariant subset of the phenomenological Ar- 
gonne i>i8 potential |22j | plus a screened Coulomb poten- 
tial. This interaction is local in coordinate space and it 
is also employed in the Green's Function Monte Carlo 
(GFMC) approach ^l|- Consequently, the use of this 
potential allows for a direct comparison of results from 
the NCSM and the GFMC methods. 

2. Non-Local in momentum space: The CDB2k inter- 
action [23 is a charge-dependent NN interaction based 
on one-boson exchange. It is described in terms of covari- 
ant Feynman amplitudes, which are non-local. Conse- 
quently, the off-shell behavior of the CD-Bonn interaction 
differs from commonly used local potentials which leads 
to larger binding energies in nuclear few-body systems. 
The newly developed N 3 LO interaction [24| is also rep- 
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resented by a finite number of Feynman diagrams. This 
interaction, however, is based on chiral perturbation the- 
ory at the fourth order, which means that it is derived 
from a Lagrangian that is consistent with the symmetries 
of QCD. It represents a novel development in the theory 
of nuclear forces. It is particularly interesting to note 
that many-body forces appear naturally already at the 
next- to-next-to- leading order (NNLO) of this low-energy 
expansion. However, in this study we use solely the NN 
part of the N 3 LO interaction. This NN potential has 
previously been applied in the NCSM approach to study 
the p-shell nuclei 6 Li and 10 B |27j . 

3. Non-Local in coordinate space: A new type of 
interaction, which respects the local behavior of tradi- 
tional NN interactions at longer ranges but exhibits a 
non-locality at shorter distances, was recently proposed 
by Doleschall et al [2f|, [2(| . The authors are exploring 
the extent to which effects of multi-nuclcon forces can 
be absorbed by non-local terms in the NN interaction. 
Their goal was to investigate if it is possible to intro- 
duce non-locality in the NN interaction so that it cor- 
rectly describes the three-nucleon bound states 3 H and 
3 He, while still reproducing NN scattering data with high 
precision. Note that all other NN interactions give a 
large underbinding of A > 3 systems. In practice, the 
INOY interaction was constructed by combining an in- 
ner (< 3 fm) phcnomenological non-local part with a 
local Yukawa tail. Hence the name INOY (Inside Non- 
local Outside Yukawa). The so called IS version of this 
interaction, introduced in Ref. |26j |. contains short-range 
non-local potentials in 1 5'o and 3 Si — 3 Di partial waves 
while higher partial waves arc taken from Argonnc vi%. In 
this study we are using the IS-M version, which includes 
non-local potentials also in the P and D waves j2j|. It 
is important to note that, for this particular version, the 
on-shell properties of the triplet P-wave interactions have 
been modified in order to improve the description of 3N 
analyzing powers. The 3 Pq interaction was adjusted to 
become less attractive, the 3 Pi became more repulsive, 
and the 3 P2 more attractive. Unfortunately, this gives a 
slightly worse fit to the Nijmegen 3 P phase shifts. 



III. RESULTS 

By construction, the ab initio NCSM method is guar- 
anteed to converge either by calculating the effective in- 
teraction using a fixed cluster approximation (e.g., two- 
body) and then solving the eigenvalue problem in increas- 
ing model spaces (A max — > oo), or by working in a limited 
model space but increasing the clustering of the effective 
interaction towards the full A-body one. Our codes are 
currently constructed to derive effective interactions up 
to the level of three-body clustering (with or without 
three-body forces). However, in this study we have cho- 
sen to approach convergence by trying to maximize our 
model space and, therefore, we limit ourselves to the use 
of two-body effective interactions. Thus we are able to 



reach the 9hQ model space for nuclei with A = 11. This 
maximal space corresponds to basis dimensions of dp = 
2.0 x 10 8 ( 9 Be), 7.0 x 10 s ( n Bc), and 1.1 x 10 9 ( n B). 
For 13 C, which is briefly discussed in connection to the 
parity-inversion problem, the largest space that we were 
able to reach was 8hQ (dp = 8.2 x 10 8 ). 

Note that model spaces with an evcn(odd) number of 
HO excitations give negative- (positive-) parity states for 
the nuclei under study. When constructing a full spec- 
trum we combine the A max M7 and (A max + 1)M7 results, 
with A max being an even number. In connection to this, 
we should also point out that very few states in 9 Be 
and 11 Be are particle stable. However, in the NCSM 
approach, all states are artificially bound due to the trun- 
cated model space and the use of HO basis functions. 

In addition to a careful study of the level ordering in 
9 Be and 11 Be, with a particular focus on the position 
of the positive-parity states, we also calculate electro- 
magnetic moments and transition strengths. For this we 
use traditional one-body transition operators with free- 
nucleon charges. Note that, due to the factorization of 
our wave function into CM and intrinsic components, we 
obtain translationally invariant matrix elements for all 
observables that we investigate in this work. However, we 
have not renormalized the operators, which means that 
the results are not corrected for the fact that we work in 
a truncated model space. The theoretical framework for 
performing this renormalization is in place, and the pro- 
cess is underway |28j . Until we are ready to implement 
the use of effective operators, we can get an indication 
on the need for renormalization by studying the basis-size 
dependence of our calculated observables . 

The cluster decomposition of the 11 Be ground state 
into 10 Bc+ n is of particular interest due to the small neu- 
tron separation energy. We have employed the formalism 
recently developed in Refs. [20l l2lj to calculate cluster 
overlap functions using our NCSM wave functions. These 
results are presented in Sec. IIII CI 



A. Dependence on HO frequency 

The first step in our study is a search for the optimal 
HO frequency. In principle, the intrinsic properties of the 
nucleus should not depend on the particular value of H£l 
since it only enters the Hamiltonian © through a CM- 
dependent term. In practice however, due to the cluster 
approximation of the effective interaction, our results will 
be sensitive to the choice of Ml. Furthermore, by con- 
struction, the effective interactions depend on the size of 
the model space, A max , and on the number of nuclcons, 
A. In order to investigate these dependences we have 
performed a large series of calculations for a sequence of 
frequencies. The results from this study are presented in 
Figs, n and [21 (for 9 Be and n Be, respectively) as curves 
showing the frequency dependence of the binding energy 
in different model spaces. We have studied this depen- 
dence for the lowest state of each parity. We are look- 
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ing for the region in which the dependence on f2 is the 
smallest; and we select this frequency (from the calcu- 
lation in the largest model space) to use in our detailed 
investigation of excited states. In our present case, this 
optimal frequency always corresponds to an energy mini- 
mum. Note, however, that the NCSM is not a variational 
method and the convergence of the binding energy with 
increasing model space is not always from above. 

Following this procedure for each nucleus and inter- 
action, we obtain the optimal HO frequencies that are 
listed in Table A few general remarks regarding the 
HO dependence, observed for the different interactions 
in Figs.Q] an d[21 can be made: (1) Clear signs of conver- 
gence is observed. The dependence on f2 becomes weaker 
with increasing size of the model space, and the relative 
difference between the calculated ground-state energies 
is in general decreasing. Furthermore, the optimal fre- 
quencies for the largest model spaces of each parity (8M7 
and 9hQ) arc approximately the same. This motivates 
our use of a single frequency to compute both positive- 
and negative-parity states; (2) This single frequency is 
found to be in the range of about Ml= 11-13 MeV for 
all interactions except for INOY, which seems to prefer a 
significantly larger HO frequency {Ml= 16-17 MeV); (3) 
In general, the behavior of the AV8', CDB2k, and N 3 LO 
interactions are very similar, but with N 3 LO having the 
largest dependence on HVt\ (4) As could be expected, since 
it is the only NN potential which is capable of reproduc- 
ing 3N binding energies, the INOY interaction exhibits 
a distinctively different behavior compared to the three 
others. The dependence on fl is encouragingly small, but 
the ground-state energy is still changing with increasing 
basis size. This is particularly true for the positive-parity 
state. We also note, from the insets of Figs. and [21 that 
it is the only interaction for which the resulting binding 
energies are approaching the experimental values. 



TABLE I: Selected optimal HO frequencies (in [MeV]). These 
choices are based on the frequency variation studies presented 
in Figs. [U and H 



Nucleus 




Interaction 






INOY 


CDB2k N 3 LO 


AV8' 


y Be 


16 


12 11 


12 


n Be 


17 


13 12 


12 



B. 9 Be 

By studying the HO frequency dependence of the 9 Be 
binding energy obtained with different NN interactions 
(see Fig. ^) it is clear that the CDB2k results have a 
slightly better convergence rate and a weaker HO fre- 
quency dependence than AV8' and N 3 LO. The INOY re- 
sults display an even weaker frequency dependence, but 
the binding energy is still moving with increasing iV max . 



It is clear from Table ITT1 that all interactions, with the 
possible exception of INOY, underbind the system. Ac- 
tually, by studying the convergence rate of the INOY 
results in Fig. ^i, it seems as if this interaction will 
eventually overbind 9 Be. This observation confirms that 
the additional binding, usually provided by 3V forces, 
can be produced by the INOY interaction. The other 
three NN interactions underbind by 12-14%. The local 
AV8' potential was also used in a recent GFMC study 
of negative-parity states in 9 Be , and we find an ex- 
cellent agreement with their ground-state binding energy 
(see Table©. 

In principle, the frequency dependence for each ex- 
cited state should be studied in order to compute its 
energy. This is particularly true in our case where we 
want to compare negative- and positive-parity states. It 
is therefore very encouraging that we find the same opti- 
mal frequency for the first positive-parity state as for the 
negative-parity ground state; and we select this frequency 
to use in our detailed investigation of excited states. In 
Figs. Hand Table HJ we present our NCSM low -energy 
spectra for different NN interactions and compare the 
results to known experimental levels. As can be seen, 
the AV8', CDB2k and N 3 LO interactions give the same 
level ordering and very similar excitation energies. It is 
noteworthy that all these high-precision NN interactions 
perform equally when applied to the A = 9 system. We 
let the AV8' spectrum shown in Fig. |3| be the graphical 
representation of all of them. Using the AV8' we can also 
make a comparison to the recent GFMC calculation [2^ . 

In general, we observe a very reasonable agreement 
with experimental levels of natural parity, while the 
unnatural-parity states are consistently high in excita- 
tion energy. For both parities, there is a general trend of 
convergence with increasing model space. When plotting 
the negative- and positive-parity spectra separately, it is 
evident that the relative level spacings are almost inde- 
pendent on the model space, so that the level ordering 
within each parity projection is remarkably stable. It is 
clear, however, that the relative position of the negative- 
versus positive-parity states is still not converged. Fur- 
thermore, when studying the AV8' convergence pattern 
in the upper panel of Fig- El it seems as if this interaction 
will predict the positive-parity states at too high exci- 
tation energies even when the calculations will be con- 
verged. This finding is consistent with an overall trend 
observed in other NCSM calculations, and it has been 
speculated whether a three-body force will correct this 
behavior H3. Although we are still not able to apply a 
true three-body force in a large enough model space, we 
get some indications from the performance of the INOY 
interaction. In Fig. 21 we see that, for this interaction, 
the positive-parity states are even higher in small model 
spaces, but that they are also dropping much faster with 
increasing N max . This issue is investigated in further de- 
tail in Sec. IIII Dl where we discuss the important question 
of parity inversion, and the general trend of the position 
of natural- versus unnatural-parity states. 



6 



-35 
-40 



> 

JB -45 



-60 



50 




• « 2M2 






■ ■ 3h£2 






•— "4h£i 






■ ■ 5hQ 


55 




• • 6h£J 






7hQ 






•— • 8MJ 






--9h£i 




'Be 



Argonne V8' ; 



10 11 12 13 



14 15 16 
h£2 



17 18 19 





FIG. 1: (Color online) The dependence on HO frequency for the calculated 9 Be^| i |j (solid lines) and 9 Be^|^~ |J (dashed 

lines) binding energies. Each panel correspond to a particular AW interaction: (a) AV8', (b) CDB2k, (c) N 3 LO, (d) INOY; 
and each separate line corresponds to a specific model space. The insets demonstrate how the minima of the curves converge 
as the model space is increased. The horizontal lines are the experimental values. 



In addition to an increase in binding energy, it has been 
found that the level ordering for many nuclei can be sen- 
sitive to the presence of multi-nucleon forces [2^, |3(], 0OJ . 
This sensitivity is the largest for those states where the 
spin-orbit interaction strength is known to play a role. 
For 9 Be we find that our calculations with the AV8', 
CDB2k and N 3 LO interactions predict the first-excited 
negative-parity state to be a 1/2 - , while experiments 
show that it is a 5/2~ (Note, however, that the CDB2k 
interaction predicts these two states to be almost degen- 
erate, and to exhibit a convergence trend indicating a 
possible level crossing at larger model spaces.). This level 
reversal was also found in the GFMC calculations using 
AV8'. The INOY interaction, on the other hand, gives 
the correct level ordering, but instead overpredicts the 
spin-orbit splitting. By performing a calculation in a 
smaller model space using the AV8' plus the Tucson- 
Melbourne TM'(99) jl^l three-nucleon interaction, wc 
found a similar result as with INOY. 

Our discussion up to this point has been concentrated 
on the low-lying levels in 9 Bc. However, in response to 
the recent evaluation published by the TUNL Nuclear- 



Data Evaluation Project |33j, we have also decided to 
summarize our results for higher excited states. It is im- 
portant to note that the experimental widths of these 
states are generally quite large, and to compute them 
correctly with the NCSM method requires a very large 
model space. Furthermore, at high excitation energies, it 
is very probable that there will be some admixture of 2MI 
intruders, and these are usually predicted too high in the 
NCSM. In any case, our results can serve as an important 
guideline as to which p-shell states that can be expected 
to appear in the spectrum, and consequently should be 
looked for in experiments. In Tabic ITTT1 wc present all 
levels that we have calculated using the CDB2k interac- 
tion in the 8M7 and 9M7 model spaces. We also show the 
tabulated experimental levels below E x = 13 MeV, taken 
from the most recent evaluation |33l | . A quick comparison 
with the previous, published evaluation [34| (from 1988), 
reveals that several new levels have been discovered and 
many spin-parity assignments have been changed. Alto- 
gether, these changes lead to a much better agreement 
with our results. In the negative-parity spectrum, our 
calculations give the correct level ordering for the first six 
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FIG. 2: (Color online) The dependence on HO frequency for the calculated u Bef |^ |J (solid lines) and n Bef |j (dashed 

lines) binding energies. Each panel correspond to a particular 7V7V interaction: (a) AV8', (b) CDB2k, (c) N 3 LO, (d) INOY; 
and each separate line corresponds to a specific model space. The insets demonstrate how the minima of the curves converge 
as the model space is increased. The horizontal lines are the experimental values. 



states. In particular, wc correctly reproduce the second 
3/2~ and 5/2~ states that were introduced in the new 
evaluation. On the other hand, we find a third 3/2 _ state 
and a second 1/2 " state that have not been observed in 
experiments. However, these states are not fully con- 
verged in our 8M7 calculation, and they are still moving 
towards higher excitation energies. We also find a 9/2~ 
state, which is fairly stable, and that has not been ex- 
perimentally identified. In the positive-parity spectrum, 
the 6.76 MeV level has now been changed to being a 
9/2+ which agrees well with our level ordering. Finally, 
it is interesting to note that our second l/2 + , 3/2 + , and 
5/2+ levels all appear in between the first 9/2 + and 7/2 + . 
None of these three states have, however, been experi- 
mentally confirmed. 

In Table UTI we also present our results for the ground- 
state quadrupole and magnetic moments, as well as 
for selected electromagnetic transition strengths. We 
should stress that free-nucleon charges have been used 
in these calculations, and that the operators have not 
been renormalized. On the other hand, the stability of 
our results can be judged by investigating the depen- 



dence on the model space. We find that the calculated 
ground-state magnetic moment and the £?(M1; | — > 

| ) transition strength are almost converged, and in 
fair agreement with the experimental values. The re- 
sults for electromagnetic quadrupole obscrvablcs arc, 
however, steadily increasing with basis size enlargement 
and should clearly benefit from the use of effective op- 
erators. As an example we can consider the evolu- 
tion of the N 3 LO results: For the {4 - 6 - S}Ml 
sequence of model spaces these observablcs increase 
as Q gs = {+3.96 - +4.10 - +4.21} [efm 2 ], and 

B (e2;§~ -> §~) = {14.9 - 15.7 - 16.7} [e 2 fm 4 ]. We 

would also like to highlight the fact that INOY gives 
much smaller values for Q gs and B(E2) than the other in- 
teractions. This is partly due to the fact that, for INOY, 
our selected frequency is much larger than for the other 
potentials which, in our limited model space, corresponds 
to a smaller radial scale. In principle, a HO frequency de- 
pendence study should be made for each operator. How- 
ever, we have also applied the INOY interaction in stud- 
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TABLE II: Experimental and calculated energies (in [MeV]) of the lowest negative- and positive-parity states in 9 Be. Quadrupole 
and magnetic moments (in [efm 2 ] and [[in]) for the ground state, as well as E2 and Ml strengths for selected transitions (in 
[e 2 fm 4 ] and [jtifv"])- Results for the AV8', CDB2k, N 3 LO and INOY AW interactions are presented. These calculations were 
performed in the 8(9)hQ model space for negative-(positive-)parity states, using the HO frequencies listed in Table Q The 
GFMC results [2t| are shown for comparison. Experimental values are from l33j[ . E x + denotes the excitation energy relative 
to the lowest positive-parity state. 



9 Be 




Exp 


INOY 


NCSM 

CDB2k 


N 3 LO 


AV8' 


GFMC 
AV8' 
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) -58.16 
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E x ( 
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1 ^ 

2 , 
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E x ( 
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j ->l 


1 ^ 

2 , 


| 2.43 


2.96 
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2.64 


2.70 


2.1 
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1 - 
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1 ' 

2 , 


| 2.78 


4.57 


2.68 


2.33 


2.50 


1.7 


E x ( 


3 - 
,2 2 


1 ' 

2 , 


| 5.59 a 


7.02 


4.98 


4.53 


4.74 




E x ( 


7 - 
s 2 i 


1 ' 

2 , 


| 6.38 


8.09 


7.80 


7.40 


7.56 


6.4 



E 



E 
l- 



i 

2 

— E e , 



E x + 
E x + 
E x + 
E x + 



1 + 

a , 

r l 

2i 2 

'3+ I 

2 , 2 , 



-56.48 
1.68 


1.37 
3.02 a 
5.08 



-50.95 
5.10 


1.39 
4.06 
6.22 



-47.81 
3.35 


1.68 
3.60 
6.36 



-47.57 
2.90 


1.68 
3.37 
6.21 



"The experimental spin-parity assignment of this level is "less 
certain" according to the TUNL Nuclear Data Evaluation I i i . 



-46.84 
3.35 


1.66 
3.49 
6.24 





5.288(38) 


3.52 


4.01 


4.21 


4.01 


5.0(3) 




-1.1778(9) 


-1.06 


-1.22 


-1.24 


-1.22 


-1.35(2) 




27.1(2) 


10.9 


14.9 


16.7 


15.0 






0.54(6) 


0.46 


0.37 


0.36 


0.37 





ies of A = 3,4 systems, for which convergence can be 
easily reached. It was found that, in particular the rms 
proton radius is always undcrprcdicted. The same result 
was obtained in Ref. j3f| through exact solutions of the 
Faddeev-Yakubovski equations, and it demonstrates that 
the interaction is too soft, resulting in a faster condensa- 
tion of nuclear matter. In this work, we have studied the 
9 Be point-nucleon radii as well as the strong El transi- 
tion from the first-excited to the ground state using the 
AV8' interaction. However, since these results will be 
compared to 11 Be data, we postpone the discussion to 

sec, mm 

In Tables II VI and [3 we show the resulting Nhfl- 
configurations and the single-particle occupancies of the 
9 Be wave function, obtained with the four different in- 
teractions. Although these quantities are not physical 
observables, they can still give interesting information. 
We see that the wave functions obtained with the AV8', 
CDB2k and N 3 LO interactions are almost identical, while 
the INOY wave function has a considerably larger frac- 
tion of low-M7 excitations. This fact is in part due to 
the higher HO frequency being used in the INOY calcu- 
lations. Furthermore, from Table llVl we see that this par- 



ticular interaction gives a different distribution of HQ ex- 
citations for the and states; which would 

indicate that the latter is slightly more deformed. How- 
ever, this behavior is not observed for the other inter- 
actions. The differences in occupations of single-particle 
levels reflect some properties of the interactions. The fact 
that the OP3/2 and 0d 5 / 2 levels have a larger occupation 
in the INOY eigenstates is direct evidence for a stronger 
spin-orbit interaction. 



C. n Be 

Most of the observations made for 9 Be concerning the 
frequency dependence of the calculated binding energies 
also hold true for 11 Be. In general, however, the sensitiv- 
ity to HQ is stronger in the n Be case. Another important 
remark, that can be made from studying Fig. [3 is that 
the binding energy of the first positive-parity state calcu- 
lated with the INOY interaction is clearly not converged. 
The relative shift in energy is actually slowly increasing 
with model-space enlargement. 
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TABLE III: Experimental and calculated energies (in [MeV]) 
of the lowest negative- and positive-parity states in 9 Be. The 
calculations were performed in the 8(9)hi} model space for 
negative-(positive-)parity states, using the CDB2k AW inter- 
action with hfl = 12 MeV. This table represents a more com- 
plete compilation of our computed levels (albeit for only one 
interaction) as compared to Table HT1 Experimental values are 
from [33ll . E x + denotes the excitation energy relative to the 
lowest positive-parity state. 



Negative parity states 
9 Be Exp CDB2k 



Positive parity states 
9 Be Exp CDB2k 



_E g , 

E.r 
E, r 
E.r 
E.r 
E. r 
E.r 
E.r 
E.r 
E.r 
E.r 
E.r 



2 

- ^ 

1- 1_ N 

2 , 2 , 



7- 1 
2, 2^ 

! 2, 
3- 1 
2 , 2 , 



- 1 
2 

- 1 
2 

- 1 
2 

5- i N 

2, 2 



-58.16 -51.16 

2.43 2.78 
2.78 2.68 
5.59" 4.98 
6.38 7.80 
7.94 a 7.96 
11.26 
11.86 
12.45 
11.28" 12.61 
11.81 13.02 




"The experimental spin-parity assignment of this level is "less 
certain" according to the TUNL Nuclear Data Evaluation 1 3 it . 
''Calculated in a smaller, model space. 
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FIG. 3: (Color online) Excitation spectrum for 9 Be calcu- 
lated using the AV8' interaction in 0hQ-9h£l model spaces 
with a fixed HO frequency of HQ = 12 MeV. The experimen- 
tal values are from Ref. |2J. The AV8' results obtained by 
the GFMC method I2H are shown for comparison (note that 
only negative-parity states were computed). The two lower 
graphs show separately the negative- and positive-parity spec- 
tra, while the upper graph shows the combined spectrum with 
selected states. 



TABLE IV: Calculated configurations of the first negative- 
and positive-parity states in 9 Be. Results obtained in our 
largest model spaces (8hQ and 9hQ, respectively) are pre- 
sented. The calculations were performed with the HO fre- 
quencies listed in Table [I] 



9Be (l 1 I) (8fin model space) 



NN interaction 


om 


2hQ 


AhQ. 


6hQ 


SKI 


INOY 


0.58 


0.19 


0.13 


0.06 


0.04 


CDB2k 


0.49 


0.22 


0.15 


0.08 


0.06 


N 3 LO 


0.47 


0.24 


0.16 


0.08 


0.06 


AV8' 


0.49 


0.22 


0.15 


0.08 


0.06 



Be (l^ I) ( gm model space) 



NN interaction 


inn 


3hQ 


5HQ, 


7Ki 


9hQ 


INOY 


0.53 


0.22 


0.14 


0.07 


0.04 


CDB2k 


0.48 


0.23 


0.15 


0.08 


0.06 


N 3 LO 


0.46 


0.24 


0.16 


0.08 


0.05 


AV8' 


0.49 


0.23 


0.15 


0.08 


0.06 



TABLE V: Calculated occupations of neutron single-particle 
levels for the first negative- and positive-parity states in 9 Be. 
Results obtained in our largest model spaces (8M1 and 9M1, 
respectively) are presented. The calculations were performed 
with the HO frequencies listed in Table 



9Be (l 1 §) ( 8m model space) 
NN interaction 0s 1/2 0p 1/2 0p 3 / 2 Ifi/g 0*4/2 0^5/2 



INOY 1.804 0.454 2.380 0.043 0.064 0.069 

CDB2k 1.773 0.511 2.277 0.067 0.060 0.072 

N 3 LO 1.768 0.521 2.256 0.077 0.060 0.072 

AV8' 1.778 0.519 2.273 0.064 0.061 0.071 



9 Be(^i+ |j (9/jfi model space) 

NN interaction 0si/ 2 Opi/2 OP3/2 lsi/2 0^3/2 0^5/2 

INOY 1.792 0.498 1.428 0.573 0.111 0.379 

CDB2k 1.761 0.530 1.377 0.666 0.113 0.317 

N 3 LO 1.757 0.536 1.365 0.676 0.117 0.312 

AV8' 1.765 0.535 1.376 0.664 0.113 0.313 
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FIG. 4: (Color online) Excitation spectrum for 9 Be calculated 
using the INOY interaction in 0hQ-9hfl model spaces with a 
fixed HO frequency of hi} = 16 MeV. The experimental values 
are from Ref. |3^l . The two lower graphs show separately the 
negative- and positive-parity spectra, while the upper graph 
shows the combined spectrum with selected states. 



FIG. 5: (Color online) Excitation spectrum for lx Be calcu- 
lated using the AV8' interaction in 0hQ-9h£l model spaces 
with a fixed HO frequency of hQ. = 12 MeV. The experimental 
values are from Ref. I39l . The two lower graphs show sepa- 
rately the negative- and positive-parity spectra, while the up- 
per graph shows the combined spectrum with selected states. 



The experimental ground state of 11 Be is an intruder 
1/2+ level, while the first p-shell state is a l/2~ situated 
at E x = 320 keV. The neutron separation energy is only 
503 keV, and there are no additional bound states. This 
level-ordering anomaly constitutes the famous parity- 
inversion problem. A number of excited states have been 
observed in different reactions and beta-decay studies. 
However, as can be seen from the summary presented in 
Table IVT1 there are considerable ambiguities in the spin- 
parity assignments. 

The low-lying experimental spectrum is compared to 
our NCSM calculated levels (obtained using the four dif- 
ferent NN interactions) in Figs. and in Table IV11I 
In the two lower panels of these figures we show, sep- 
arately, the negative- and positive-parity spectra, while 
the upper panel shows the combined spectrum with se- 
lected states. Note that, those experimental levels for 
which there is an uncertainty in the parity assignment, 
are included in all three panels. There are clear signs 
of convergence with increasing model space. However, 
as was also observed for 9 Be, the relative position of the 
negative- and positive-parity spectra has clearly not con- 
verged, and the latter is still moving down. The most 



dramatic drop is observed in the INOY spectrum, thus 
indicating the importance of a 3iV force. With this par- 
ticular interaction, the l/2 + level actually ends up below 
all, but one, of the negative-parity states in the (8-9) Ml 
calculation. We refer to Sec. lIIIDl for further discussions 
on this topic. 

We stress again that the relative level spacings, ob- 
served when plotting negative- and positive-parity states 
separately, is remarkably stable. Furthermore, the or- 
dering of the first six(four) levels of negative (positive) 
parity, is the same for all four potentials. This calcu- 
lated level ordering is summarized in Table IVIIII Our 
results can, therefore, provide input to help resolve the 
uncertainties of the experimental spin-parity assignments 
(cf. Tabic rvT|l . Note in particular that some experiments 
suggest that there are three low-lying 3/2~ states. The 
task to compute three levels with the same spin quickly 
becomes very time consuming with increasing dimension, 
since it requires many Lanczos iterations. Therefore, this 
third state was studied in two separate runs, using only 
the AV8' and INOY potentials, and is included in Figs.0 
and |H1 up to the QMl model space. These calculations 
confirm the existence of three low-lying 3/2~ levels, but 
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TABLE VI: Present situation of the spin-parity assignments for the lowest states in Be. The table contains published results 
from the FAS evaluation of 1990 |3fil | and from more recent experimental studies. These studies include direct reactions such 
as (t,p) (Liu-Fortune) and 12 C( 11 Be, 11 Be / ) (Fukuda) in which the extracted angular distributions were analyzed using DWBA 
theory. The remaining references are measurements of /3-delayed neutrons in coincidence with 7-rays. All decays that were 
observed in these experiments had \og(ft) values that were consistent with allowed transitions, indicating that the corresponding 
final states have negative parity and J < §. 
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TABLE VII: Experimental and calculated energies (in [MeV]) 
of the lowest negative- and positive-parity states in 11 Be, as 
well as the magnetic moment (in [^tjv]) of the ground state. 
Results for the AV8', CDB2k, N 3 LO and INOY NN inter- 
actions are presented. These calculations were performed in 
the 8(9)ftf2 model space for negative-(positive-)parity states, 
using the HO frequencies listed in Table Q E x - denotes the 
excitation energy relative to the lowest negative-parity state. 
Experimental values are from l39l. l42| . 



n Be 




NCSM 




Exp 


INOY CDB2k N 3 LO AV8' 



E 



(21 2) 



Err 



E x - 
E x - 
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0.32 


? a 
? a 
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-2.44 -2.54 





2.08 
3.70 
4.43 





2.04 
3.81 
4.38 




-65.48 -59.51 

? a 1.85 



-54.26 


2.01 



-54.13 -52.98 

1.98 2.03 



/'gs 



-1.6816(8) -1.47 -1.55 -1.58 -1.58 



"There are large ambiguities in the experimental spin-parity as- 
signments, cf. Table IVTI 



TABLE VIII: NCSM observed ordering (from left to right) of 
11 Be negative- and positive-parity states (separately). Note 
that all four NN interactions used in this study give the 
same ordering for the first six(four) negative-(positive-)parity 
states. 
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FIG. 6: (Color online) Excitation spectrum for lx Be calcu- 
lated using the CDB2k interaction in 0h£l-9hQ model spaces 
with a fixed HO frequ ency of hQ. = 13 MeV. The experimental 
values are from Ref. 139^1 . The two lower graphs show sepa- 
rately the negative- and positive-parity spectra, while the up- 
per graph shows the combined spectrum with selected states. 



they also stress the presence of a 5/2~ state which is not 
completely consistent with Rcfs. [3(| H3, HI] . However, 
we can not rule out the possibility of a low-lying intruder 
2fi,0-dominated state, which would avoid detection in our 
study. These states have a different convergence pat- 
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FIG. 7: (Color online) Excitation spectrum for 1 Be calcu- 
lated using the N 3 LO interaction in 0fiQ-9hQ model spaces 
with a fixed HO frequency of HQ = 12 MeV. The experimental 
values are from Ref. I39L The two lower graphs show sepa- 
rately the negative- and positive-parity spectra, while the up- 
per graph shows the combined spectrum with selected states. 



FIG. 8: (Color online) Excitation spectrum for lx Be calcu- 
lated using the INOY interaction in 0hQ~9hi} model spaces 
with a fixed HO frequ ency of ftfi = 17 MeV. The experimental 
values are from Ref. |39l . The two lower graphs show sepa- 
rately the negative- and positive-parity spectra, while the up- 
per graph shows the combined spectrum with selected states. 



tern than Ohfl states and generally appear at too high 
an excitation energy in the smaller model spaces, sec e.g. 
Ref. [13. 

In summary, our results suggest that there are two 
excited positive-parity states below 4 MeV (rather than 
three as stated in Ref. [Hi). The 1.78 MeV level should 
be a 5/2+ state, while either the 3.41 or the 3.89 MeV 
level is a 3/2+ . Our results do not support the pres- 
ence of a high-spin (J > 7/2) state, which one can find 
in Ref. [3|j. We do observe three low- lying 3/2~ states 
although they arc accompanied by a 5/2 - state. 

The strength of the electric dipole transition between 
the two bound states in 11 Be is of fundamental impor- 
tance. This is an observable which has attracted much 
attention since it was first measured in 1971 |40j . and 
again in 1983 HJ. The cited value of 0.36 W.u. is still 
the strongest known transition between low- lying states, 
and it has been attributed to the halo character of the 
bound-state wave functions. Unfortunately, by working 
in a HO basis, we suffer from an incorrect description 
of the long-range asymptotics, and we would need an 
extremely large number of basis states in order to re- 
produce the correct form. This shortcoming of the HO 



basis is illustrated by the fact that we obtain a value 
for the El strength which is 20 times too small (see Ta- 
bic [Kj . When studying the dependence of this value on 
the size of the model space, we observe an almost linear 
increase, indicating that our result is far from converged. 
For the {(4-5) - (6 - 7) - (8 -9)} Ml sequence of 

model spaces, the 11 Bc El strength, B ^El; |^ — > f^J) 

calculated with the AV8' interaction increases as: 
{0.0054 - 0.0059 - 0.0065} [e 2 fm 2 ]. The correspond- 
ing sequence of results for 9 Be is: B ^El; ^ — ► ^ 

{0.029 - 0.031 - 0.033} [e 2 fm 2 ], which demonstrates a 
similar increase. However, for this nucleus we note that, 
in the largest model space, our calculated El strength is 
only off by a factor of two compared to experiment . In 
addition, a consistent result is found for the much weaker 
| + — > I El transition in 9 Be, where we also obtain a 
factor of two smaller 73(E1) than experiment. These re- 
sults accentuates the anomalous strength observed for 
11 Be. A simple explanation for the failure of HO calcula- 
tions in the 11 Be case was given by Millcncr et al fllj. It 
was shown that there is a strong cancellation in the cal- 
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culated El transition amplitude due to the insufficient 
description of the long-range asymptotics (see in partic- 
ular Tables IV and V in Ref. [41|). By simply replacing 
their HO single-particle wave functions with solutions to 
the Schrodinger equation with a Woods-Saxon potential, 
they found that the magnitude of the neutron ls!/ 2 0pi/2 
single-particle matrix element increased significantly so 
that the cancellation was removed. Even though our 
multi- Ml calculations give a significant improvement of 
the calculated El strengths as compared to their simple 
(0—l)Ml model, the underlying problem is still present. 

Another operator which is sensitive to the long-range 
behavior of the wave function is the point-nucleon ra- 
dius. However, even though no operator renormalization 
has been applied, our results show a fair stability with in- 
creasing model space, and they are in rather good agree- 
ment with experimental findings for both 9 Be and 11 Be 
(see Table IIX|I . It is probably safe to assume that the 
missing part of the 11 Be matter radius originates mainly 
in an underestimation of the point-neutron radius. One 
should also remember that the experimental results for 
matter radii, in these light systems, are highly model- 
dependent and are usually theoretically extracted from 
measurements of the interaction cross section. In addi- 
tion, we have also calculated the radii of the first excited 
state. For both isotopes it is found that the unnatural- 
parity state has a 10% larger neutron radius than the 
natural-parity one, probably due to a larger admixture 
of sd-shell neutrons. Finally, the ground-state magnetic 
moment of 11 Be has been measured |42j and we find a 
reasonable agreement with our calculated value, see Ta- 
ble eiii 

TABLE IX: Nuclear ground-state radii (in [fm]) and the El 
strengths (in [e 2 fm 2 ]) for the strong ground-state transitions 
in 9 Be and n Be. The NCSM calculations were performed in 
the 8(9)ftfi model space for negative-(positive-)parity states 
using the AV8' interaction. The GFMC result for 9 Be, with 
the same interaction |29[|. is shown for comparison. Experi- 
mental values are from |33l l39t l49t l5d . 





Rn 


Rp 


9 Be( 

^mat 


B(E1) 

1+ 3- 5+ 3- 


Exp 
NCSM 
GFMC 


2.40 


2.39 
2.27 
2.41(1) 


2.45(1)° 
2.34 


0.061(25) 0.0100(84) 
0.033 0.0057 




Rn 


Rp 


n Be( 

^mat 


i+ s\ 

B(E1) 

1- 1 + 

2, 2, 


Exp 
NCSM 


2.66 


2.30 


2.86(4) 
2.54 


0.116(12) 
0.0065 



"Interaction radius 



The standard halo picture of the 11 Be ground state 
is a simple two-body configuration consisting of an inert 
10 Be core coupled to an S1/2 valence neutron. Theoret- 



ical estimates of the spectroscopic factor for this com- 
ponent range from 0.55 to 0.92, see e.g. Table 1 in 
Ref. 0] • The experimental situation is also unclear since 
the extracted results are generally model-dependent. In 
the literature one can find values from 0.36 to 0.8, see 
e.g. Fig. 8 in Ref. j^. An important question is to 
which extent the first-excited 10 Be(2+) state contributes 
to the simple two-body configuration. The formalism 
for investigating cluster structures of NCSM eigenstates 
was recently developed in Ref. [2l]]. We have calcu- 
lated the overlap of the n Bc^i^ |^ state with different 

10 Be + n channels. To this aim, the 11 Be( 10 Be) wave 
functions were calculated using the CDB2k interaction 
in a 7(6)M) model space. We used a HO frequency of 
Ml = 14 McV, which corresponds to the optimal value for 
calculating binding energies in these two model spaces. 
The largest overlap functions (in jj coupling) are pre- 
sented in Fig. while the corresponding spectroscopic 
factors (the overlap function squared and integrated over 
all r) arc summarized in Table ^ Several additional 
channels, such as the overlap with the second excited 
2^ state in 10 Be, were also computed but their spectro- 
scopic factors were found to be very small (< 0.001). 
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FIG. 9: (Color online) The largest radial overlap functions for 

the n Be^i + |J state decomposed as 10 Be + n in jj-coupling. 

The results presented here were obtained with the CDB2k 
interaction (Ml = 14 MeV) with the 11 Be( 10 Be) wave function 
calculated in a 7(6)M1 model space. The thin, dotted line 
shows the dominant overlap function calculated in a smaller 
5(4)ftfi model space. 

Several observations can be made when studying these re- 
sults: (1) The 11 Be ground state has a large overlap with 
[ 10 Be(0^) <g> n (si/ 2 )] (5" = 0.82), but also with the core- 
excited [ 10 Be(2^) ® n (d 5/2 )] channel (S = 0.26). These 
results are in good agreement with the consensus of re- 
cent experimental studies, see e.g. Refs. 0,0. (2) The 
thin dotted line in Fig.0shows the [ 10 Be(0^) <g> n (s 1/2 )] 
overlap function calculated in a smaller model space. 
From this comparison it is clear that the results are quite 
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TABLE X: Spectroscopic factors for the n Be |j ground 

state decomposed as 10 Be + n in coupling. The results 
presented here were obtained with the CDB2k interaction 
(HQ = 14 MeV) with the n Be( 10 Be) wave function calcu- 
lated in a 7(6)hQ model space. For comparison, we list spec- 
troscopic factors extracted from three recent experiments uti- 
lizing different reactions. 
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Be ® n 




Transfer 


Knockout 


Breakup 
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{1,3) 


NCSM 


[43j a 


[45J 6 


[44] c 




(o.j) 


0.818 
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0.78 


0.61(5), 0.77(4) 


2i 


2,| 


0.263 


0.09-0.16 










0.022 










0.032 










0.005 








o 8 + 


o,| 


0.037 









"DWBA analysis of "Befed). 
6 From 9 Be ( 11 Be, 10 Bc + 7) 

Spectroscopic factors extracted from 11 Be breakup on lead and 
carbon targets respectively. 



stable with regards to a change in iV max . The interior 
part does hardly change at all, while the tail is slowly ex- 
tending towards larger inter-cluster distances. This state- 
ment is true for all channels shown in the figure except 
for those involving the two high-lying + states [see bul- 
let (4) below]. (3) The inset shows the main component 
plotted on a logarithmic scale. This graph clearly demon- 
strates the fact that our HO basis is not large enough to 
reproduce the correct asymptotic behavior. Even though 
the tail is extending further with increasing W max , it still 
does not reach the expected exponential decay. Instead 
it dies of too fast. (4) Our calculated 10 Be O4" and 0^ 
states are found to be 2MI dominated, and their binding 
energies have not converged in the NCSM calculation. 
The cluster overlaps with these states do not display the 
same stability as observed for the other channels. In- 
stead, there is a large dependence on W max . A similar 
result was found in Ref. [21j and it is just another man- 
ifestation of the slower convergence of the 2MI states in 
the NCSM. 

Finally, we compare, in Tables IXll and IxTTl the result- 
ing configurations and the occupancies of single-particles 
states obtained with different interactions. Again, it is 
clear that the INOY eigenstates have a larger fraction of 
\ow-Ml excitations and that this interaction results in a 
different single-particle spectrum due, most likely, to a 
stronger spin-orbit interaction. 

D. Parity inversion 

One of the main objectives of this study has been to 
investigate the relative position of negative- and positive- 
parity states in the region around 11 Bc. As we have 
shown, none of our calculations reproduce the parity in- 
version that is observed for this nucleus. However, con- 



TABLE XI: Calculated configurations of the first negative- 
and positive-parity states in 11 Be. Results obtained in our 
largest model spaces {8T1Q. and 9hQ, respectively) are pre- 
sented. The calculations were performed with the HO fre- 
quencies listed in Table 



llBe (^ I) (8ft^ model space) 
AW interaction OhQ. 2HD. AhQ, 6hQ. 8HQ, 



INOY 0.59 0.17 0.14 0.06 0.04 

CDB2k 0.51 0.20 0.15 0.08 0.06 

N 3 LO 0.49 0.22 0.15 0.08 0.06 

AV8' 0.48 0.21 0.16 0.08 0.07 



n Be(^i+ f J (9/ift model space) 
WW interaction lhQ 3htt 5fi,fi 7hQ 9HQ 
INOY 056 O20 014 O06 O04 

CDB2k 0.50 0.21 0.15 0.08 0.06 

N 3 LO 0.49 0.22 0.16 0.08 0.06 

AV8' 0.48 0.22 0.16 0.08 0.07 



TABLE XII: Calculated occupations of neutron single-particle 
levels for the first negative- and positive-parity states in 11 Be. 
Results obtained in our largest model spaces (8M7 and 9hQ, 
respectively) are presented. The calculations were performed 
with the HO frequencies listed in Table 

n Be(±~ (8HZI model space) 
WW interaction 0s 1/2 0pi/ 2 OP3/2 lsi/2 0^3/2 0^5/2 



INOY 1.862 1.078 3.643 0.046 0.065 0.075 

CDB2k 1.835 1.093 3.597 0.066 0.062 0.072 

N 3 LO 1.832 1.095 3.586 0.073 0.062 0.072 

AV8' 1.828 1.094 3.579 0.073 0.061 0.072 



n Bc(i+ |j (9hQ model space) 

WW interaction 0s 1/2 Opi/2 OP3/2 lsi/2 0^3/2 0^5/2 

INOY 1.845 0.504 3.300 0.658 0.086 0.345 

CDB2k 1.824 0.600 3.181 0.742 0.088 0.285 

N 3 LO 1.823 0.616 3.153 0.752 0.091 0.281 

AV8' 1.820 0.630 3.135 0.768 0.088 0.265 



sidcring the slower convergence rate for lfiil-dominated 
states in the NCSM, and the large, but still finite, model 
spaces that we were able to use, our results are actu- 
ally very promising. In all nuclei, we found a fast drop 
of the unnatural-parity states with respect to the natu- 
ral ones. This behavior has already been demonstrated 
in earlier NCSM studies, but the drop that we observe 
in 11 Be is the most dramatic so far. Furthermore, the 
results obtained with the INOY interaction are clearly 
different from the others, which indicates the signifi- 
cance that a realistic 3N force should have in a funda- 
mental explanation of the parity inversion. Note that 
INOY is a two-body interaction, but that it simulates 
the main effects of 3A^ forces by short-range, non-local 
terms. Furthermore, the 3 P NN interactions are slightly 
modified in order to improve the description of 3N an- 
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alyzing powers. Fig. ^| shows the calculated excitation 
energy of the first positive-parity states in 9 Be and 11 Be 
as a function of the basis size, iV max . For illustrative pur- 
poses we have extrapolated our results to larger model 
spaces assuming an exponential dependence on V max , 
i.e., E x = E XtOC + aexp (— bN max ). Note that the (0- 
1)M7 points are excluded from the fits. The extrapolated 
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FIG. 10: (Color online) Basis size dependence of the cal- 
culated E x (| + ) excitation energy relative to the lowest 

negative-parity state in (a) 9 Be and (b) 11 Be. The results 
for four different NN interactions are compared. For each po- 
tential, a single, fixed HO frequency was used (see Table 0. 
The dashed lines correspond to exponential fits of the calcu- 
lated data and, for illustration, these curves are extrapolated 
to larger model spaces. 

INOY results end up below the other interactions; and 
for 9 Be the curve is actually approaching the experimen- 
tal value. With all other interactions, the extrapolated 
excitation energy is wl-2 MeV too high. 

When discussing the position of the first unnatural- 
parity state, it is very interesting to study the system- 
atics within the A = 11 isobar and the N = 7 iso- 
tonc. To this aim, we have performed large-basis cal- 
culations for 11 B and 13 C. The diagonalization of the 
11 B Hamiltonian in the 9M7 space proved to be our 
largest calculation so far. For 13 C we were only able 
to reach the 8M7 space. Both studies were performed us- 



ing the CDB2k interaction and an HO frequency of Ml = 
13 MeV. The ground-state binding energies (obtained in 

the 8hn space) are: E ( n B; §~ §) = 66.25 MeV, and 

E ( 13 C; \~ ±) = 86.53 MeV. Our calculated U B spec- 
trum, including the first negative-parity state for each 
spin up to J = 9/2, plus the lowest positive-parity state, 
is compared to known experimental levels in Fig. 1111 
Note that we obtain an incorrect 1/2" ground-state spin 
in our largest model space. However, the first 3/2~ and 
1/2" states are found to be almost degenerate, and there 
is a trend indicating that the position of these levels may 
eventually intersect as the basis size is increased. In prin- 
ciple, a thorough frequency variation study should be 
performed in order to clarify the fine details of the pre- 
dictions. In any case, it is clear that the level splitting is 
described incorrectly with this interaction. Basically the 
same result was found in an earlier NCSM study |31j us- 
ing a three-body effective interaction derived from AV8'. 
In that paper, it was also shown that the correct level 
ordering can be reproduced, and the splitting greatly im- 
proved, by adding a realistic 3N force. For 13 C we have 
only computed the lowest state for each parity. However, 
this nucleus has also been studied previously using the 
NCSM. A spectrum obtained with the CDB2k interac- 
tion was presented in Ref. |4^, while calculations with 
a genuine 3N force were reported in Ref. [3l||. In both 
papers, the study was limited to negative-parity states 
and a smaller model space {AMI) was used. 
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FIG. 11: (Color online) Excitation spectrum for 1 H calcu- 
lated using the CDB2k interaction in 0ftf2-9M7 model spaces 
with a fixed HO frequency of hQ = 13 MeV. The experimental 
levels are from Ref. [39ll . Note that there are many additional 
levels between the experimental l/2 + and 9/2" shown in the 
figure. However, we have only computed the first level for a 
given spin, and the l/2 + was the only positive-parity state 
that was considered. 

Let us now comment on our n B and 13 C results and 
return to the important question of the position of the 
first positive-parity state. The calculated l/2 + excita- 
tion energy for these two nuclei, as a function of V max , is 
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shown in Fig. ^] It is a fascinating empirical fact that, 
by simply going from Z = 4 — > Z = 6, the first 1/2+ 
state moves from being the ground state in 11 Be to be- 
come an excited state at 3.1 MeV in 13 C. In the odd-Z 
nucleus n B, the first positive-parity state is found quite 
high in the excitation spectrum, namely at 6.8 MeV. It 
is a significant success of the NCSM method, and of the 
iViV interactions being employed, that these huge shifts 
are accurately reproduced in our calculations. However, 
as can be seen from Figs. 1101 and !12l the calculated exci- 
tation energy always turns out to be too large. A com- 
parison of our extrapolated CDB2k results shows that 
they exceed the experimental values by ssl-2 MeV for 
all four isotopes. As a final remark, our INOY results for 
9 Be and n Bc indicate that the use of a realistic 3V force 
in a large basis space might correct this discrepancy 
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FIG. 12: (Color online) Basis size dependence of the cal- 
culated E x excitation energy relative to the lowest 

negative-parity state in (a) n B and (b) 13 C. These results 
are obtained with the CDB2k interaction and hQ — 13 MeV. 
The dashed lines correspond to exponential fits of the calcu- 
lated data and, for illustration, these curves are extrapolated 
to larger model spaces. 



IV. CONCLUSIONS 

We have performed large-basis ab initio no-core shell 
model calculations for 9 Be and 11 Be using four different 
realistic AW interactions. One of these, the non-local 
INOY interaction, has never before been used in nuclear 
structure calculations. Although it is formally a two- 
body potential, it reproduces not only AW data (beware 
of the fact that the 3 P interactions are slightly modified 
in the IS-M version that we are using) but also the bind- 
ing energies of 3 H and 3 Hc. Therefore it has been of 
particular interest for our current application, where we 
have striven to maximize the model space by limiting our- 
selves to AW interactions, but have still been very much 
interested in the effects of three-body forces. We have 
computed: binding energies, excited states of both par- 
ities, electromagnetic moments and transition strengths, 
point-nucleon radii, and also the core-plus-neutron clus- 
ter decomposition of the n Be ground state. 

In summary, for the calculated spectra we found clear 
signs of convergence, and a remarkable agreement be- 
tween the predictions of different AW interactions. In 
particular, the relative level spacings observed when plot- 
ting positive- and negative-parity states separately, were 
found to be very stable and to agree well with experimen- 
tal spectra. This has allowed us to make some conclusions 
regarding the largely unknown spin-parities of unbound, 
excited states in n Be. An overall observation is that the 
AV8' and N 3 LO potentials produce very similar results, 
while CDB2k gives slightly more binding. The INOY in- 
teraction is clearly different; giving a much larger binding 
energy and a stronger spin-orbit splitting. Both these ef- 
fects would be expected from a true 3N force, but are 
here achieved by the use of short-range, non-local terms 
in the AW interaction. 

Furthermore, it was also clear from our study that 
our results for observables connected to long-range op- 
erators, have not converged. These calculations would 
clearly benefit from operator renormalization, in order 
to correct for the limited model space being used. In 
particular, the extremely strong El transition between 
the two bound states in 11 Be, was underestimated by a 
factor of 20. We have discussed how this illustrates the 
fact that the anomalous strength is due to the halo char- 
acter, and hence large overlap, of the initial and final 
state wave functions; a property which is extremely hard 
to reproduce using a HO basis. In the NCSM approach, 
there is no fitting to single-particle properties, e.g., by 
the use of empirical interactions. Instead, the effective in- 
teractions are derived from the underlying intcr-nuclcon 
forces. Therefore, it is likely that a good description of 
loosely bound, and unbound, single-particle states might 
require a very large number of HO basis functions. 

An important topic of this work has been the investi- 
gation of the parity inversion found in 11 Be. We did not 
reproduce the anomalous l/2 + ground state in our ab 
initio approach, but did observe a dramatic drop of the 
positive-parity excitation energies with increasing model 
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space. Furthermore, the behavior of our INOY results 
suggests that a realistic 37V force will have an important 
influence on the parity inversion. However, in order to 
pursue this question further, an improved computational 
capacity is needed. We have also performed large-basis 
calculations for 11 B and 13 C. In this way, we were able to 
put our 11 Be positive-parity results into a wider context 
by studying the systcmatics within the Z = 4 isotopes 
( 9 Be), the N = 7 isotonc ( 13 C), and the A = 11 isobar 
( n B). Although we found that the NCSM always over- 
estimates the excitation energy of the first unnatural- 
parity state, we did reproduce the very large shifts ob- 
served for these different nuclei. This is an important 
finding which leads us to the optimistic conclusion that 
the parity-inversion problem should be possible to repro- 



duce in the NCSM starting from realistic inter-nuclcon 
interactions. 
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